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■ Abstract. We discuss the case of massive gravitons and their relation with the 

cosmological constant, considered as an eigenvalue of a Sturm-Liouville problem. A 
^ ' variational approach with Gaussian trial wave functionals is used as a method to 

i study such a problem. We approximate the equation to one loop in a Schwarzschild 

background and a zeta function regularization is involved to handle with divergences. 
The regularization is closely related to the subtraction procedure appearing in the 
computation of Casimir energy in a curved background. A renormalization procedure 
is introduced to remove the infinities together with a renormalization group equation. 
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1. Introduction 

There are two interesting and fundamental questions of Einstein gravity which have not 
received an answer yet: one of these is the cosmological constant Ac and the other one 
H ' is the existence of gravitons with or without mass. While the mass-less graviton is a 
natural consequence of the linearized Einstein field equations, the massive case is more 
delicate. At the linearized level, we are forced to introduce the Pauli-Fierz mass term^] 

= if / ^'^V^-^ [^''^ V - h'] , (1) 

where is the graviton mass and k = SnG. G is the Newton constant. The Pauli-Fierz 
mass term breaks the symmetry h^^, — > h^y -|- 2V (^^S,u)^ but does not introduce ghosts. 
Boulware and Deser tried to include a mass in the general framework and not simply in 
the linearized theory. They discovered that the theory is unstable and produce ghostsj21. 
Another problem appearing when one consider a massive graviton in Minkowski space is 
the limit mg — *• 0: the analytic expression in the massive and in the mass-less limit does 
not coincide. This is known as van Dam-Veltman-Zakharov (vDVZ) discontinuity 
Other than the appearance of a discontinuity in the mass-less limit, they showed that 
a comparison with experiment, led the graviton to be rigorously mass- less. Actually, 
we know that there exist bounds on the graviton rest mass that put the upper limit on 
a value less than 10~^^ — 10~^^(7|[4j. Recently there has been a considerable interest in 
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massive gravity theories, especially about the vDVZ discontinuity examined in de Sitter 
and Anti-de Sitter space. Indeed in a series of papers, it has been shown that the vDVZ 
discontinuity disappears in the mass-lees, at least at the tree level approximation [Hj, 
while it reappears at one loopjB]. If we fix our attention on the positive cosmological 
term expanded to one loop, we can see that its structure is 



Ac 



(2) 



which is not of the Pauli-Fierz form|. Nevertheless, we have to note that the non trace 
terms of Sp,f. and Sa^ can be equal if 

= A,. (3) 



In other words the graviton mass and the cosmological constant seem to be two aspects 
of the same problem. Furthermore, the cosmological constant suffers the same problem 
of smallness, because the more recent estimates on A^ give an order of 10~'^'^ GeV^ , while 
a crude estimate of the Zero Point Energy (ZPE) of some field of mass m with a cutoff at 
the Planck scale gives Ezpe ~ lO^^GeT^^ with a difference of about 118 orders jH]. One 
interesting way to relate the cosmological constant to the ZPE is given by the Einstein 
field equations withour matter fields 



R 



■fiu 



(4) 



where G^i, is the Einstein tensor. If we introduce a time-like unit vector such that 



u ■ u 



-1, then 



(5) 



This is simply the Hamiltonian constraint written in terms of equation of motion. 
However, we would like to compute not Ac, but its expectation value (Ac) on some 
trial wave functional. On the other hand 



-R + 
2k 



where R is the scalar curvature in three dimensions. Therefore 



(Ac) 



1 

V 



—A / d^xK 



1 

'v 



— ( / d-^xA^ 



where the last expression stands for 

lJV[g,,]^* [g,,]J^d'xn^[g,, 



l(^ 


f^d^xAj, 




V 







A 



(6) 



(7) 



(8) 



V JV[g,^]^*[g,,]^[g,,] 

and where we have integrated over the hypersurface E, divided by its volume and 
functionally integrated over quantum fluctuation. Note that Eq.dH)) can be derived 
starting with the Wheeler-De Witt equation (WDW) jH] which represents invariance 
under time reparametrization. Extracting the TT tensor contribution from Eq.(|H)) 

I To this purpose, see also Ref.[7]. 
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approximated to second order in perturbation of the spatial part of the metric into 
a background term, gij, and a perturbation, hij, we get 



The propagator K-*- (x, a;)^^^; can be represented as 



2A fr) 



(9) 



(10) 



where h[l^'^ (it) are the eigenfunctions of A2. r denotes a complete set of indices and 
A (r) are a set of variational parameters to be determined by the minimization of Eq.®. 
The expectation value of is easily obtained by inserting the form of the propagator 
into Eq.Q and minimizing with respect to the variational function Aj (r). Thus the 
total one loop energy density for TT tensors is 



A(Ai) 



4^ 



ul (r) + Jul (r) 



The above expression makes sense only for (r) > 0. To further proceed, we count the 
number of modes with frequency less than Ui, i = 1,2. This is given approximately by 



Ui{l,uOi) (2/ + 1) 



(12) 



where Vi{l^uji), i = 1,2 is the number of nodes in the mode with {l,uji), such that 
(r = r (x)) 



1 f+°° I 

Pi (/, ^i) = — / dxJ kf (r, /, uJi). 
27r y_oo ^ 



(13) 



Here it is understood that the integration with respect to x and I is taken over those 
values which satisfy kf {r,l,Ui) > 0, i = 1,2. Thus the one loop total energy for TT 
tensors becomes 

2 



1 



E 

1=1 



dx 



9 l^i 
duj 



-duJi 



(14) 



2. The massive graviton transverse traceless (TT) spin 2 operator for the 
Schwarzschild metric and the W.K.B. approximation 

The further step is the evaluation of Eq. fll4|) . when the graviton has a rest mass. 
Following R,ubakov|12j. the Pauli-Fierz term can be rewritten in such a way to explicitly 
violate Lorentz symmetry, but to preserve the three-dimensional Euclidean symmetry. 
In Minkowski space it takes the form 

Sm= / d'^Xy/^Cm, (15) 
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where 

= mlh^^hoo + 2mlh^%oi - mjh'^hij + mlh''hjj - 2mlh^^hii (16) 

A comparison between Sm and the Pauh-Fierz term shows that they can be set equal if 
we make the following choice§ 

2n 22222^n 
TTIq = U 171^ = 1712 — ""^3 — ""^4 = ^ > U. (17) 

If we fix the attention on the very special case = mf = m"^ = m\ = 0; = > 0, 
we can see that the trace part disappears and we get 

Its contribution to the Spin-two operator for the Schwarzschild metric will be 

(A,/.^^)^ := - [ArK'^)] + 2 {Rh^^X + i^yni (19) 

and 

- = -A, (ft-)f + ^ (l - Hi^) (h--)] . (20) 

A5 is the scalar curved Laplacian, whose form is 

^ _/ 2MG\ (f ^/2r-3MG\ d 
\ r ) dr"^ \ J dr 

and Rj is the mixed Ricci tensor whose components are: 

I lyiJ ijniJ ly'^J I 

This implies that the scalar curvature is traceless. We are therefore led to study the 
following eigenvalue equation 

(As/i^^)^' = u^h] (23) 

where uj"^ is the eigenvalue of the corresponding equation. In doing so, we follow 
Regge and Wheeler in analyzing the equation as modes of definite frequency, angular 
momentum and parity ^UJ- In particular, our choice for the three-dimensional 
gravitational perturbation is represented by its even-parity form 

{h^''^-)] (r, ^, 0) = diag [H (r) , K (r) , L (r)] Y,^ (^, 0) . (24) 

Defining reduced fields and passing to the proper geodesic distance from the throat of 
the bridge, the system becomes 

(25) 

-£ + ^-^ + ml{r)\Mx)=u;lMx) 
See also Dubovski|13| for a detailed discussion about the different choices of toi, m2, 771,3 ^md 771,4 
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where we have defined r = r (x) and 

ml (r) = + [/i (r) = mg + mj (r, M) — 7713 (r, M) 

(26) 

m| (r) = m'^g + U2 (r) = m'^g + (r, M) + m| (r, M) 

(r, M) ^ when r — » 00 or r — » 2MG and (r, M) = 3MG/r^. Note that, while 
m2 (r) is constant in sign, mf (r) is not. Indeed, for the critical value f = 5MG/2, 
ml (f) = and in the range {2MG,5MG/2) for some values of m^, (f) can be 
negative. It is interesting therefore concentrate in this range, where ml (r, M) vanishes 
when compared with m^ (r, M). So, in a first approximation we can write 

mj (r) ~ — m2 (ro, M) = — mg (M) 

(27) 

m2 (r) ^ + m2 (ro, M) = + (M) 

where we have defined a parameter ro > 2MG and mg (M) = SMG/rg. The main 
reason for introducing a new parameter resides in the fiuctuation of the horizon that 
forbids any kind of approach. It is now possible to explicitly evaluate Eq. fjl4|) in terms 
of the effective mass. One gets 

K I 

A = Pi + P2 = -— ^ 2I / . ^V""'' ~ ^''^'^'^ (2^) 

where we have included an additional 47r coming from the angular integration. 



3. One loop energy Regularization and Renormalization 



Here, we use the zeta function regularization method to compute the energy densities p\ 
and p2- Note that this procedure is completely equivalent to the subtraction procedure 
of the Casimir energy computation where the zero point energy (ZPE) in different 
backgrounds with the same asymptotic properties is involved. To this purpose, we 
introduce the additional mass parameter [i in order to restore the correct dimension for 
the regularized quantities. Such an arbitrary mass scale emerges unavoidably in any 
regularization scheme. Then we have 

1 



Pi 



2£ 



duji- 



{ijj} — m'i 



The integration has to be meant in the range where 



m: 



I (r) > 0. One gets 



Pi {£) 



K- 



mf (r) 

256772 



1 , 
- + In 

e 



P' 



mf (r) 



+ 21n2- i 
2 



(29) 



(30) 



i = 1,2. To handle with the divergent energy density we extract the divergent part of 
A, in the limit e — and we set 

G 



A 



div 



32716 



(m^ (r) + m2 (r)) 



(31) 
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Thus, the renormahzation is performed via the absorption of the divergent part into the 
re-definition of the bare classical constant A 



A ^ Ao + A"^™. 



(32) 



The remaining finite value for the cosmological constant reads 



Ar 



1 



SttG 2567r2 



+ni2 (r) 



In 



m\ (r) 



In 



\m\ (r)| 



+ 2 In 2 - 



1 



m\ (r) 



+ 21n2 - - 
2 



(Pi (P) + P2 (p)) = Peff (P> r) .(33) 



The quantity in Eq. ipHjl depends on the arbitrary mass scale /x. It is appropriate to 
use the renormahzation group equation to eliminate such a dependence. To this aim, 
we impose that[TT] 



1 (/i) _ d_ 



P-rPeffiP^r). 



(34) 



Solving it we find that the renormalized constant Aq should be treated as a running one 
in the sense that it varies provided that the scale /i is changing 



G 



P 



Ao (/i, r) = Ao {jJo, r) + —- (m^ (r) + (r)) In 

lOTT Ho 



Substituting Eq.§^ into Eq.(^ we find 

Ao{po,r) 1 



SttG 



2567r2 



+ {ml + ml{M)y 



In 



{ml - ml (M)) 
ml + ml (M) 



In 



m^ — ml (M) I 

pI 



(35) 



21n2 + - 
2 



pI 



2 In 2 + 



(36) 



We can now discuss three cases: 1) m^ 3> ml (M), 2) m^ = ml (M), 3) <^ ml (M) 
In case 1), we can rearrange Eq. fl36|) to obtain 



Ao {po,r) 



m„ 



In 



WiJ 2 



SvrG 1287r2 
where we have introduced an intermediate scale defined by 



pIi = pI exp 



3mg (M) 
2m4 



(37) 



(38) 



With the help of Eq. ()38|) . the computation of the minimum of Aq is more simple. Indeed, 
if we define 

,2 ^.,4 



m„ 



X 



In ix) + 



1 



(39) 



^Pm ^ 
As a function of x, Kq^m{pqiX) vanishes for x = and x = exp ( — |) and when 
X G [O, exp (— I)] , Ao,M {Po, x) > 0. It has a maximum for 



1 

X = - 
e 



m 



2 _ ^p\i _ ^Po 



exp 



3m^ (M) 
2m4 



(40) 
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and its value is 



Ao,M (yWo, X) = 7^3^ = ^—^ exp 



or 



27re^ 



G 



3m^ (M) 



Ao,M (/^o, a;) = exp 

In case 2), Eq-dHUj) becomes 

Ao(/io,r) Ao(/io) 



3m^ (M) 



or 



SttG SttG 
Ao (/io) (M) 



1287r2 



In 



SvrG 



In 



1" 



+ 



+ 



Again we define a dimensionless variable 



m„ 



X 



Ao,o (/"o, x) 



Gil 



In (x) + - 



(41) 



(42) 



(43) 



(44) 



(45) 



4/i^ SttG tt 

The formal expression of Eq. ()45|1 is very close to Eq. ()39|1 and indeed the extrema are in 
the same position of the scale variable x, even if the meaning of the scale is here different. 
Ao,o (/^O) x) vanishes for x = and x = 4 exp (— |) • In this range, Ao,o (/^o; a;) > and it 
has a minimum located in 



and 



or 



1 

X = - 
e 



Ao,o (/io,a;) 



Ao,o (/Uo,a;) 
Finally the case 3 ) leads to 



27re2 

^rn\ = ^rnt (M) . 
327r 3 327r °^ ^ 



Ao (/io, r) (M) 



SttG 



1287r2 



In 



(M) 
4/i^ 



+ 



where we have introduced another intermediate scale 



/^m = /^o exp 



3m„ 



2m^ (M) 

By repeating the same procedure of previous cases, we define 



X 



Ao,m (/io, x) 



4 



TT 



Also this case has a maximum for 
1 



4/y2 4//2 



In fx) 



3m 



(46) 

(47) 
(48) 

(49) 

(50) 

(51) 
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and 

A iu -X)- - exD ( (53) 

or 

G ( "irn^ \ 

Ao,M (yt^o, = — (M) exp . (54) 

Remark Note that in any case, the maximum of A corresponds to the minimum 
of the energy density. 

A quite curious thing comes on the estimate on the "square graviton mass" , which 
in this context is closely related to the cosmological constant. Indeed, from Eq. ()46|) 
applied on the square mass, we get 

oc /i^ ~ lO^^GeV^ = lO^^eV^, (55) 

while the experimental upper bound is of the order 

(m^^^^oclO-^«-10-^«e\/^ (56) 

which gives a difference of about 10^^ — 10^°^ orders. This discrepancy strongly recall the 
difference of the cosmological constant estimated at the Planck scale with that measured 
in the space where we live. 
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